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ABSTRACT 

Wc have used the European VLBI Network (EVN) to observe a sample of Lyman 
j^ ■ Break Analogs (LB As), nearby (z < 0.3) galaxies with properties similar to the more 

distant Lyman Break Galaxies (LBGs). The study of LBGs may help define the feed- 
back relationship between black holes (BHs) and their host galaxies. Previous VLA 
observations have shown that the kpc-scale radio emission from LB As is dominated by 
starbursts. The main targets of this VLBI experiment were selected because they pos- 
sessed emission-line properties between starbursts and Type 2 (obscured) AGN. Eight 
^^ ■ targets (three star-forming LBAs, four composite LBAs, and one Type 1 AGN) were 

tH- I observed at 5 GHz, four of which (one star-forming LBA and three composite LBAs) 

were also observed at 1.7 GHz. One star- forming LBA and one composite LBA were 
detected above 5cr at 1.7 GHz (only), while the AGN was detected at 5 GHz. In both 
LBAs, the radio luminosity (Ln) exceeds that expected from supernovae (remnants) 
based on a comparison with Arp220, Arp229A and Mrk273, by factors of 2 — 8. The 
composite LBA exhibits a compact core emitting around 10% of the VLA flux density. 
The high Tb of 3.5 x 10^ K and excess core L^ with respect to the Ln/Lx relation of 
radio-quiet AGN indicate that this LBA possesses an obscured AGN {Mbh ~ 10^"^ 
Mq). In three other composite LBAs detected previously in the X-ray, no radio sources 
d i were detected, indicating either variability or the presence of an obscured AGN below 

our radio sensitivity. While weak AGN may co-exist with the starbursts as shown in 
at least one of the LBAs, their contribution to the total radio flux is fairly minimal. 
Our results show that the detection of such weak AGN presents a challenge at radio. 
X-ray and optical emission-line wavelengths at z ^ 0.2, indicating the great difficulties 
that need to be overcome in order to study similar processes at high redshift when 
these types of galaxies were common. 

Key words: radio continuum: galaxies, techniques: interfcromctric, galaxies: active, 
galaxies: starburst, galaxies: ISM 



1 INTRODUCTION lowing the mass increase of their host galax ies in accordance 

• ,,,,, with the Mbh-o or Mt,h-M t,uiae relations (|Magorrian et al.1 

1 he nuclei oi local galaxies contain supermassive black holes. L^onl 1^ i i — ,. — '. — nr;7UU]\ n . i r n i i , 

^ , ^, , ^^^ , ^ ., 119981 : IGebhardt et al.ll2000l l. Furthermore, feedback between 

Evidence suggests that these BHs have grown over time rol- , , „tt i . i • i . • , ■ i i , r , 

these BHs and their hosts is an essential element ot structure 
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formation models as it has been proposed that they regulate 
or suppress star formation as a result of active galactic nu- 
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cleus (AGN) activity fe.g. lDi Matteo et al.ll2005h . One must 
thus approach galaxy formation with both star formation 
and BH accretion in mind. 

Starburst galaxies at high redshift, such as the UV- 
selected Lyman Break Galaxies or their dusty counterparts 
such as the sub-mm galaxies (SMGs), are the precursors of 
present day massive galaxies undergoing a p hase of intense 
star formation (e.g. lAdelberger et al.lbOOa '). A significant 
fraction of their stellar mass is contained in large, dense stel- 
lar clumps. It has been suggested that these giant clumps 
could host seed BHs th at merge in the center of the galaxy 
(|Elmegreen et al.ll200a ). where they can continue to grow 
further through other accretion events. It has been shown 
that the majority of SMGs at high redshift host an AGN a s 
indicated by their X-ray emission ({Alexander et al.ll2005h . 
while the fraction a mong LBGs is a few percent at most 
llSteidel et all |2002|: iLehmer et all l2005l : iLaird et all l2006l : 
lOuchi et al.ll2008l : IZheng et al.ll201Cl ) . One possible explana- 



tion for this discrepancy is a delay expected between the mo- 
ment of onset of the starburst and the moment that a galaxy 
is able to feed its BH, eis one has to wait for the strong winds 
and supernovae (SN) from the most m assive stars to have 
ceased before gas can fa ll onto the BH (|Norman fc Scovilld 
ll988l : lDavies et al.| l2007''). This would be consistent with the 
relatively young ages of typical LBGs at z ~ 3 — 4 compared 
to the dustier and more evolved SMGs. 

Although the majority of LBGs do not appear to host 
luminous AGN, it is difficult to detect the presence of a 
low-luminosity AGN (LLAGN), especially when it is at the 
center of a powerful starburst. From studies in the local uni- 
verse we know that weak or absent X-ray emission does not 
necessarily imply that an accreting BH is absent. LLAGN 
have a large range in optical-to- X-ray luminosity ratios, and 
man y are very weak in har d X-rays due to heavy obscura- 
tion (JHeckman et al.ll2005al ') . Other AGN identification tech- 
niques, such as the commonly u sed diagnostic tests b ased on 
optical emission line ratios (e.g. [Baldwin et al.lll98ll . "BPT" 
diagrams) can also yield inconclusive results. These line ra- 
tios are susceptible to, e.g., dust , meta l licity, age and sta r 
formation history (|Kewlev et al.l l2006l : lYuan et all |2010| ). 
New low velocity shock models furthermore show that the 
line ratios arising in purely star-forming or starbursting en- 
vironments are often indistinguishable from those that are 
due to a combination of star formation and some AGN con- 
tribution, both oc cupying the sam e "composite" region of 
the BPT diagram (IRich et aLlbOllI ). These problems can be 
partly alleviated by disentangling the emission from differ- 
ent galaxy components (e.g., nucleus vs. outer parts) using 
integral field spectroscopy, but this is currently not feasible 
at high redshifts except for the most extended galaxies. 

We previously selected a rare population of nearby 
{z < 0.3) galaxies that are remarkably similar to LBGs 
in most of their basic physical properties. These "Lyman 
Break Analogs" (LBAs) are similar to LBGs in mass, age, 
size, metallicit y, optical extinction, SFR, and gas veloc- 



ity dispersion ( Heckman et al.ll2005bl : iHoopes et al.ll2007h . 
They also have morpholog ies similar to those at z = 2 — 4 



llOverzier et al.ll200 
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have similar emission line kine- 



matics (JBasu-Zvch et al.ll2009l : iGoncalves et al.ll2010f ). and 



any evidence for an unobscured (Type I) AGN (largely by 
selection), it is unclear whether some may host obscured 
(Type II) LLAGN. A small subset of LBAs falls in the com- 
posite region of the BPT diagram. These sources were also 
found to have the densest and m ost massive nuc l ei rem inis- 
cent of extremely young bulges (jOverzier et al.ll2009l '). An 
obvious question to ask is whether we see any evidence for 
BH growth. 

At low intrinsic radio luminosities, the total radio emis- 
sion of AGN is often dominated by the star formation in 
their host galaxies (e.g. IPadovani et al.ll201ll : lKimball et al.l 
I2OIH ). High-resolution Very Long Baseline Interferometry 
(VLBI) observations can be used to search for the elusive 
LLAGN through their compact, high brightness tempera- 
ture radio emission amidst other compact sources corre- 
sponding to radio supernovae (RSN) and supernova rem- 
nants (SNR) frequently seen on pc scales by allowing us to 
identify excess radi o emission above that expected for RSN 
and SNR mode ls ("Kcwlcv et al] I2OO0I: iBondi e t al.' '2005|: 
Gallimore et al.1 ,2006: Pcrcz -Torres et aLlboioT TBiggs ct ^ 



a similar interstellar medium dominated b y strong feed- 
back from massive, clumpy starburst regions ( Overzier et al.l 
I2OO9I : iHeckman et al.ll2011h . While none of the sources show 



20101 ) . Previous VLBI surveys have primarily targeted dusty 
galaxies (e.g. Luminous or Ultra-Lu minous IR galaxies , 
(U)LIRGs) in the nearby universe (e.g. ICorbett et al.ll2003l : 
IParra et al.ll2010^ . The UV-selected LB A sample allows us 
to study the relationship between nuclear starbursts and 
BHs in an environment much more similar to that expected 
at high redshift. We present the results from a pilot project 
using the European VLBI Network (EVN) exploring for the 
first time the compact radio structures in the nuclei of LBAs. 
The structure of this Paper is as follows. In §2 we present 
our sample and give details on the 6 and 18 cm radio obser- 
vations and data reduction methods. In §3 we present the 
main results. We discuss the implications of our results in 
§4, and conclude with a summary in §5. Throughout this pa- 
per, we assume a cosmology [TIa{,Qa,Ho] ~ [0.27,0.73,73.0] 
(with Ho in km s^^ Mpc^^) so that the angular scales at 
z—0.1, 0.2, and 0.3 are f»1.8, 3.2, and 4.3 pc nias^^. 



2 SAMPLE, OBSERVATIONS, AND DATA 
REDUCTION 

For this experiment we targeted LBAs that were previously 
detected at the >0.5 mjy level in the Very Large Array 
(VLA) observations of lBasu-Zvch et all (|2007h . Our VLBI 
sub-sample consists of three sources that are unlikely to host 
an AGN, and four sources that may be of starburst-AGN 
composite type based on their optical emission line ratios 
(jOverzier et al.ll2009l ). For comparison, we also include an 
eighth source showing strong evidence for pure AGN emis- 
sion. We note that this source has similar far-UV properties 
as the LBAs but it does not satisfy the full LBA selection 
criteria because the latter were designed to exclude pure 
AGN. Details on the sources are given in Table [T] 

Two separate sets of observations were completed us- 
ing the European VLBI Network (EVN) over the course of 
two years. EVN observations of J0150-I-1308, J0213+1259, 
J0921-I-4509 and J2103-0728 were carried out at 1.7 GHz on 
three days between 4-8 June 2008. The participating tele- 
scopes were the Lovell Telescope in Jodrell Bank, the 14- 
element Westerbork Synthesis Radio Telescope, Effelsberg, 
Onsala (25m), Medicina, Noto, Torun, Urumqi, Shanghai 
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and Cambridge. Baselines therefore ranged from 198 km to 
8476 km. The aggregate bitrate per telescope was 1024 Mbps 
using 2-bit sampling and 4s integration time. There were 
8x16 MHz subbands each with 32 channels in both LCP 
and RCP polarizations. 

To extend the coherence time, the targets were phase- 
referenced to nearby calibrators with separations of 0.5-2 
degrees selected from the VLBA Calibrator Lisllj. The cal- 
ibrator positions were accurate to within 0.20 to 1.32 mas. 
Within the phase-reference cycle of about eight minutes, 100 
seconds was spent on calibrators and 280 seconds on the tar- 
gets, and there were 30s gaps included for system tempera- 
ture measurements. The total on-source time per target was 
about 3.5 hours, that would ideally lead to a thermal noise 
below 10 /ijy/beam. 

Follow-up EVN observations of an extended sam- 
ple, now including J0808-f3948, J0815-h5004, J1029-H4829, 
J1353-I-6648 were carried out at 5 GHz on 19-20 March 
2010 (see Table [T|. A higher frequency was used in this sec- 
ond run because it was believed to lead to higher detection 
rates: the expected flat spectrum emission from an AGN core 
would dominate over the steep spectrum emission from star- 
forming regions and it was expected that free-free absorp- 
tion would not play a dominant role at this frequency. The 
array included the Lovell Telescope, Yebes (40m), Wester- 
bork, Effelsberg, Onsala (25m), Medicina, Noto and Torun. 
The resulting baselines ranged from 266 to 2279 km. The 
data were recorded at 1024 Mbps in a similar fashion to the 
first observations (2-bit sampling and a 2-second integration 
time). There were 8x16 MHz subbands each with 32 chan- 
nels in both LCP and RCP polarizations. Targets were again 
phase-referenced to calibrators with angular separations of 
0.5-2 degrees. A phase-reference cycle of 90-300 seconds on 
the calibrators and targets respectively was employed. The 
total on-source time on the targets was between 1-2 hours. 
The expected rms noise was below ~15 /xjy/beam allowing 
for > 5cr detections for targets at the level of 100 /xJy. 

The data were correlated by the EVN Data Proces- 
sor at the Joint Institute for VLBI in Europe (JIVE). For 
data reduction and analysis we used the NRAO Astronom- 
ical Image Processing System (AIP S). We followed stan- 
dard data reduction procedures (e.g. iDiamond Ill995l ): the 
data were amplitude calibrated with measured gaincurves 
and system temperatures then fringe-fitted and bandpass 
calibrated using the compact calibrator sources. Anoma- 
lous solutions were manually flagged in AIPS using the task 
snedt to view calculated rate and delay solutions. Subband 
8 for the Lovell telescope and subband 7 for Torun were 
completely flagged in the second set of observations us- 
ing task uvflg due to poor bandpass solutions. The cal- 
ibrated target data were averaged in time and frequency 
and then exported to UVFITS for mat. Additional manua l 
flagging was carried out in Difmap (jShepherd et al. Ill994h . 
For the second set of data collected on March 2010, ampli- 
tude calibration was verifled on several compact calibration 
sources (J0209-H352, J0808-F4950, J1344-F6606) in Difmap. 
The sources were cleaned using repeated calls to clean and 
selfcal and then the amplitude corrections were examined 
using gscale. It was noted that across all calibrator sources 
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all Torun amplitudes were off by a factor of 1.7 while Noto 
subband 7 amplitudes were off by a factor of 1.5 and sub- 
band 8 amplitudes were off by a factor of 1.7. Systematic 
corrections were made to these subbands using the AIPS 
task clcor and then the data was re-exported to UVFITS 
format. Possible sources of such systematic errors include 
faulty Tsys measurements due to RFI usually present across 
most subbands. 

Dirty maps were formed with natural weighting (uvw 0, 
-2) using the mapplot command in Difmap. For each target 
source the peak brightness and rms noise was determined. 
The results are listed in Tabled For the two safe detections 
we show clean maps in Figure [l] 



3 RESULTS 

For the 1.7 GHz data two sources, J0150H-1308 (star- 
forming LBA) and J092H-4509 (composite LBA) were de- 
tected above 5a. J0150-fl308 is located at {a,5)j2ooo of 
(01:50:28.39531,-^13:08:58.0833). J092H-4509 is located at 
(09:21:59.40425,-1-45:09:12.4232). Both have a positional er- 
ror of approximately 1 mas. For the 5.0 GHz data only one 
source was detected above 5a, J1029+4829 (the AGN) while 
the two sources detected at 1.7 GHz were not. J1029-f4829 
is located at (10:29:59.95930,-1-48:29:38.0012) with a posi- 
tional error of approximately 2 mas. Real detections should 
exhibit the synthesized beam pattern in their images due to 
the beam pattern's large effect on the noise in any image, 
although this cannot be seen at very low signal-to-noise ra- 
tios. The two LBAs detected at 1.7 GHz exhibit this beam 
pattern in their dirty images suggesting they are safe detec- 
tions. The strong detection of these two sources at 1.7 GHz 
demonstrates the success of our phase-referencing efforts. 
We believe that the sources that remained undetected at 
1.7 GHz are diffuse radio sources which lack a compact core 
detectable on the scale of our observations. Applying a uv- 
taper to downweight longer baselines, testing the effects of 
possibly over-resolving the flux, did not lead to a detection 
thus strengthening our reasoning that these are indeed faint 
sources. 

At 5 GHz one likely detection just above 5 a was made. 
For the source J 1029-1-4829 (AGN) the peak brightness lo- 
cation did not vary when applying a strong UV-taper or 
when removing baselines or observation times suggesting a 
probable detection. Observations with a greater sensitivity 
would be needed to confirm this detection. The two sources 
detected at 1.7 GHz were not detected at 5 GHz. This result 
could be explained in several possible ways. One option is 
that our phase-referencing did not work well. Alternatively, 
the sources did not have compact emission well exceeding 
the lOO^Jy level so the sources were completely resolved in 
this observation. Lastly, it is possible that they had spectra 
of (a ^ —0.1, S oc v") calculated using the 5 a upper limit. 

For the detected sources, parameters were obtained by 
fitting point and/or circular-Gaussian model components to 
the uv-data using the model fit command in Difmap. The 
results are summarized in Table [S] It should be noted that 
in all cases less than 50% of the measured VLA flux density 
was recovered. For the source J0150-(-1308 (star- forming) 
only an extended model significantly larger than the restor- 
ing beam produced a good fit. Fitting a point source model 
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Table 1. Basic source parameters of our sample. The c oordinates are from S PSS. The 1.4 GHz VLA flux de nsities and rad i o SFR s are 
from Basu-Zych et al. (2007). Opt ical line ra i ^ios ar e from lOverzier et al.l | |2009| ). Lpiji and SFR{UV+IR) from lOverzier et al.l boill ). and 
2-10 kcV X-ray luminosities from ljia et al.l 1 120111 ). The radio-derived SFRs agree very well with those based on the UV-I-IR, indicating 
that the radio emission is dominated by star formation. 
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Figure 1. Naturally weighted clean maps of J0150+1308 (left) and J0921+4509 (right) at 1.7 GHz. The data for J0150+1308 were 
additionally weighted with a Gaussian taper of 0.5 at a uii-distance of 10 MA. The contour levels are ±1, 1.5, 2, 2.5 and 3 times the 
3-cr noise of 45 /ijy/beam, the tapered beam is 41.7x27.5 mas at PA 16.8 degrees, the peak brightness is 148 /ijy/beam. In case of 
J0921+4509 the contour levels arc ±1, 1.5, and 2 times the 3-cr noise of 45 /ijy/beam, the beam is 15.9x5.7 mas at PA 16.3 degrees, the 
peak brightness is 111 /-ijy/beam. 



on the map left residual emission on the dirty map, indicat- 
ing a more extended structure. Thus the non-detection of 
J0150+1308 (star-forming) at 5 GHz is most likely due to 
over-resolution. On the other hand, J0921+4509 (compos- 
ite) could be fit equally well with a point source, a relatively 
compact circular Gaussian (~ L32 mas), or a more extended 
Gaussian source model (with size about the same as the ma- 
jor axis of the restoring beam). 

Brightness temperatures were determined for each 
source detected using the formula 



Tb.vibi = 1.22 X 10''(l + 2) 



S 



Oidiu'^ ' 



(1) 



where S is the flux density (Jy), 9i and 02 are the major and 
minor axes of the fitted Gaussian model component (mas or 
FWHM for a point source) and ly is the observing frequency 
(GHz). 

Because we could not distinguish between point and 
extended models in the case of J0921+4509 (composite), we 
used a Monte Carlo simulation of 1000 trials to attempt to 
constrain the source size. All Monte Carlo simulation s were 
carried out using ParselTongue (jKettenis et al. 120061 ) to call 
the relevant AIPS task repeatedly. Taking the 1.7 GHz field 
of J0213+1259 (composite) as a non-detection field we used 
the AIPS task uvsub to add a point source directly to the uv- 
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Table 2. Measured value of rms noise, peak brightness and signal to noise ratio of the peaks for our eight targets (four targets were 
observed both at 1.658 and 4.990 GHz). Upper limits are given for non-detections. Note that because both observations used a different 
array configuration and the L-band run included Chinese and South African telescopes the resolution from the two runs was quite similar. 



Source 


V 


noise 


peak 


S/N 


beam 


Tb 




(GHz) 


(MJy/b) 


(/^Jy/b) 




(mas) 


(lO^K) 


J0150+1308 


1.658 


15 


149 


9.9 


8.6x6.2 


>32.0 


J0213+1259 


1.658 


14 


<54 


<3.9 


5.9x7.8 


<5.0 


J0921+4509 


1.658 


15 


113 


7.5 


5.7x15.9 


>1.1 


J2103-0728 


1.658 


17 


<73 


<4.3 


6.3x10.2 


<5.4 


J0150+1308 


4.990 


26 


<111 


<4.3 


6.1x8.5 


<1.1 


J0213+1259 


4.990 


26 


<115 


<4.4 


6.1x8.7 


<1.2 


J0808+3948 


4.990 


35 


<162 


<4.6 


5.2x10.1 


<1.7 


J0815-I-5004 


4.990 


31 


<117 


<3.8 


5.1x9.3 


<1.1 


J0921+4509 


4.990 


25 


<100 


<4.0 


5.4x8.1 


<1.1 


J1029+4829 


4.990 


24 


127 


5.3 


7.1x9.4 


>1.2 


J1353+6648 


4.990 


25 


<98 


<3.9 


5.5x8.3 


<1.0 


J2103-0728 


4.990 


27 


<111 


<4.1 


7.1x10.2 


<0.8 



Table 3. Modelfit results and brightness temperatures for all de- 
tections near or above 5(7. J0150+1308 (star-forming) could be 
fit well only by a circular-Gaussian model but J0921-I-4509 (com- 
posite) gave similarly good results when fit by a point source, a 
compact circular-Gaussian smaller than the beam, and a more ex- 
tended model. Here we give the compact circular-Gaussian model 
size (see text). For the weakest detection, J1029-I-4829 (AGN), 
only a point model was used. In this case the geometric average 
of the beam major and minor axes is given as upper limit for the 
size. 



Source 



J0150-I-1308 
J0921-I-4509 
J1029-I-4829 



S 



Ratio 



(GHz) (/ijy) (EVN/VLA) (mas) (XQ^VC) 



1.658 
1.658 
4.990 



456 
111 
126 



0.304 
0.079 

N/A 



59.3 
1.32 
<9.4 



0.10 

34.9 

>0.12 



data. Random noise was simulated by placing the source at 
a different location in the field with every call to uvsub. The 
AIPS task omfit was then used to fit a circular Gaussian 
model to this point source. Trials were conducted with initial 
modelfitting parameters of both 2 mas and 10 mas for the 
major axis and giving the correct location of the simulated 
source. In both cases, the results were strongly peaked at 
a size of mas though an initial fit parameter of 10 mas 
occasionally led omfit to give results between 10 mas and 
15 mas as the estimation of the size of our point source. We 
repeated the trials adding phase errors of 10 rms by using 
ParselTongue to manually edit the SN table of our dataset. 
When fitting a 10 mas model to this data the secondary peak 
between 10 mas and 15 mas was more pronounced. So while 
one of our modelfit results for J0921-f4509 had a radius of 
up to 15.9 mas our simulations reveal the radio source could 
still be unresolved with a much smaller radius, especially 
given expected phase errors in our data. 

Next we performed a similar trial but attempted to fit 
a Gaussian model to an inserted Gaussian source of 10 mas. 
The initial model provided to omfit exactly matched the 
simulated source. While results of the modelfitting program 
tended to cluster around 10 mas, there was a spread of 
±5 mas with and without phase errors of 10 rms. Thus any 



modelfit result obtained is likely accurate to within ±5 mas. 
However, this trial never produced results in the region of 
1 mas. In our simulations an extended model comparable to 
the beamsize can be fit by a weak point source accurately, 
but it is not possible to fit well an extended source with 
a compact model much smaller than the beam. Since the 
VLBI structure of J0921+4509 (composite) can be fit well 
with compact models, it is most likely that the true size of 
the source is about 1.32 mas, with an error of ±5 mas. 

To summarize, we conclude that we have made two safe 
detections, one of which, J0150-f 1308 (star- forming), was re- 
solved at 59.3 mas while the other, J0921+4509 (composite), 
is practically unresolved at 1.32 mas. In the next section we 
will discuss the interpretation of these results. 



4 DISCUSSION 

4.1 Preexisting evidence for hidden AGN in 
composite LBAs 

In order to address the question of whether weak AGN 
are present in the nuclear starburst regions of the LBAs, 
we specifically targeted four LBAs having optical emission- 
line spectra intermediate between pure starbursts and Type 
II AGN (J0213+1259, J0808-f3948, J092H-4509, J2103- 
0728). We also targeted three objects that, at least accord- 
ing to their main optical emission line ratios, are not sus- 
pected of having an obscured AGN present (J0150+1308, 
J0815-1-5004, J1353-I-6648), and one source that is an AGN 
(J1029+4829). These sources are shown in the optical emis- 
sion line diagnostics diagram in Figure [2] (BPT diagram). 

Previous work has shown that the LBAs that lie in the 
composite part of the BPT diagram have a number of other 
properties as well that make them peculiar with r espect to 
the "non-composite" LBAs: lOverzier et al.l (|2009l ) showed 
that they are host to particularly dense and compact star- 
bursts (radii of order 100 pc), reaching stellar mass densities 
comparable to those of dense star clusters, while being up to 
1000 times more massive (stellar masses of ~ 10* — lO'' Mq). 
They also tend to be the most massive galaxies in the LBA 
sample, and have the highest SFR (several tens to a hundred 
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Figur e 2. Optical emissio n lines diagnostic diagram accord- 
ing to [Baldwin et all l|l98ll . BPT). The LEA sample studied in 
lOverzier et alj 1 120091 ) is indicated by red/blue squares. Filled, 
blue squares mark LBAs that have been observed as part of our 
VLBI campaign, while (open) red squares mark LBAs that were 
not observed with VLBL Objects detected during our EVN ex- 
periment are marked by the large, open blue squares. The black 
filled square indicates the AGN-like object J1029+4829 that was 
observed as part of this sample but is not an LBA. S mall dots in- 
dicate LIRGs from the Northern COLA sample from lParra et alJ 
bOlOl ) also observed with VLBI, with detections marked by the 
large open circles. Objects that lie in the region of the diagram in 
between the dotted and dashed lines are so-called starburst-AGN 
composite objects, while objects to the left of the dashed line and 
to the right of the solid line a rc Hll-li kc and AGN-like objects, 
respectively IjKauffmann et al.li2003i : iKewlev et al.ll2006l ) . 



M(7) y r~\ the exception being J0808-l-3948'). lHeckman et all 
(J201ll ) presented FUV spectra obtained with HST/COS that 
showed that these dense, compact starburst regions drive 
powerful winds with median and maximum velocities several 
times higher than that seen in typical starbursts, including 
LBGs at high redshift. Based on the starburst nature of the 
FUV spectra, the high wind speeds and the lack of evidence 
of unobscured (i.e. Type I) AGN, it was argued that these 
are momentum-driven winds caused by the intense super- 
nova activity in the compact starbursts at the base of the 
wind. Finally, Ijia et al.l (|2011f ) found the first plausible ev- 
idence for the presence of low luminosity, obscured AGN 
in some or all of 6 composite LBAs observed in the X-ray. 
The presence of an AGN was inferred from elevated levels 
(factors of 2-10) of hard X-ray emission over that expected 
from the SFR-Lx sequence formed by pure star-forming 
galaxies (see Figure |3]). However, the SFR-Lx correlation 
and its scatter are not particularly well constrained, mak- 
ing the detection somewhat ambiguous. A stacked spectrum 
showed the potential detection of the 6.4 keV iron emission 
line (which is a key feature of obscured AGN), but only at 
the ~2(T level. |jia et al.l l|2011f ) further found that the ra- 
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Figure 3. Hard X-ray lum inosity vs. far infrared luminosity of 
LBAs from I Jia et al.l (|201l|), who observed six composite LBAs 
with XMM-Newton. The four composite LBAs also observed 
with VLBI are shown in blue, while the two composite LBAs 
not observed with VLBI are marked red. The LBA composite 
J0921-I-4509 that was detected during our VLBI experiment is 
marked with a large blue square. Other samples s hown for com- 
pariso n are as follows: LIRGs and LIRGs AGN from lLehmer et al.l 
1I2OIOI ) (open and filled triangles, resp.), LIRGs A GN (diamonds) 
and H ard X-ray Quiet galaxies (small points) from llwasawa et al.l 
1I2OO9I I. The s olid line indicates th e empirical Lx,2—l0k£V ~ ^FIR 
relation from lRanalli et al.l l|2003h . with the dashed lines showing 
an associated uncertainty of a factor of two. The dash-d otted line 
is the best-fit relation for the composite LBAs found by I Jia et al.l 
l|201lh . 



tios of the mid-infrared (24 fj,m) continuum to [O IIIJA5007 
luminosities in these LBAs was higher than the values for 
typical Type 2 AGN by an average of 0.8 dex. Combining 
all these clues, it is likely that a hidden AGN is present in 
some or all of these composite LBAs, but that their bolo- 
metric luminosities are completely dominated by the intense 
starburst. 



4.2 New insights from the VLBI data 

The VLBI data on LBAs presented in this paper provide 
another piece of information on the issue of AGN in LBAs. 
As shown in Figure [J] we have detected one out of the four 
composite LBAs observed with VLBI, and one out of the 
three Hll-like LBAs. (Our third detection of a known AGN 
will not be discussed further here). 

On first sight, these results do not offer a clear di- 
chotomy between star-forming and composite LBAs. Be- 
cause all LBAs have extended radio structures due to star 
formation and supernova activ ity extending to larger scales 
than those probed by VLBI (|Basu-Zvch et al.l [JOOtI ) . the 
emission we have detected could in principle be coming from 
certain compact, high surface brightness regions related to 
those structures. The fact that three of the four composite 
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Table 4. Radio luminosity, physical size and resulting star forma- 
tion rate for each of the 5 a detections. Should the radio emission 
be from a starburst the radio luminosity can be used to estimate 
the star formation rate. We calculated the SF R of our sources 
using the equation from iBasu-Zvch et al.l 1120071) . 



Source 


Ll.7 

(1029 ergs s-i Hz-i) 


Physical Size 
(pc) 


SFR 

(Moyr-l) 


J0150+1308 
J0921+4509 


2.116 
1.388 


146.2 

4.76 


7.78 
5.10 


J1029+4829 


1.540 


<33.7 


5.66 



LBAs do not show any VLBI radio emission at all indicates 
that (1) the radio emission arising from their star-forming 
regions is completely resolved at milliarcsecond or parsec 
scales, and (2) no (bright) radio AGN is present in these 
sources. 

What other constraints do we have in order to study 
the relative contributions from starburst and AGN activ- 
ity in this sample? Our detections have brightness tempera- 
tures of ~ 10^ and ~ 10^ K (see Table [3|. These values are 
equal to (in the case of J0150-I-1308) and well in excess of 
(J0921-I-4509) the typical upper-limit for thermal emission 
arising from star- forming regions (~ 10"" K), thus indicating 
that the emission in these objects is likely of non-thermal 
origin , at least in the case of J092H-4509 (jCondon et al.l 
Il99lh . Even though at these frequencies starburst emission 
may be dominated by nonthermal emission as well, their 
maximum brightness temperature still does not appear to 
exceed that limit. One interpretation is that the observed 
radio emission is associated with radio supernovae (RSNe) 
or supernova remnants (SNRs). The non-detection at 5 GHz 
of both of our detections at 1.7 GHz supports this conclu- 
sion: our sources have steep spectra which may be indicative 
of RSNe or SNR. 

We can compare the total (non-thermal) radio luminos- 
ity detected from J0150-fl308, L1.7 ghz = 2.1 x 10^^ erg s"^ 
Hz-^ and J0921-F4509, L1.7 ghz = 1.4 x 10^'' erg s'^ Hz"^ 
(Table |4]), with that expected based on nearby starburst 
galaxies for which it is possible to determine the relative 
contributions of individual RSNe and SNRs. The brightest 
RSN observed in starburst galaxies have a typical luminos- 
ity at 18 cm of the order of 2 x 10^^ erg s~^ Hz~^, while the 
most luminou s RSN ever observ ed (SN1988z) had 2 x 10^* 
erg s~^ Hz~^ (jBondi et aI.ll2005l V It is thus clear that the ra- 
dio emission from LBAs requires large complexes of such SN 
(or an AGN). We will employ the well-studied examples of 
such complexes observed in the s tarburst g alaxies Arp299- 
AJPcrcz- Torres ct al. 200d l2010l ). Ar p220 jLonsdale et all 
l200g r iParra et al.lbOOTi '). and Mrk273 (JBondi et al.ll2005l ') to 
guide our discussion. 

The central 150 x 85 pc region of Arp299-A consists of 
a cluster of 26 RSNe and SNRs detected at 5 GHz, with 
a total luminosity coming from this region of L5.0 ghz = 
2x10^* ergs"^ Hz~\ corresponding to Li. 7 ghz ~ 4.5x10^** 
erg s~^ Hz~^ assuming a median spectral index of -0.72 
for the region as a whole. The Western nucleus of Arp220 
is another example of such a prolific supernova factory: it 
consists of 18 individual RNSe/SNRs within a region of 50 x 
150 pc and a total radio luminosity of I/1.7 ghz ~ 1.6 x 10^® 



erg s~^ Hz~^|f| In order to minimize potential systematic 
biases when we compare these objects with the LBAs, we 
will apply a scaling based on the ratios of the physical areas 
{A) observed and the total IR luminosities (Lm) of these 
systems. The expected 18 cm luminosities for the LBAs are 

then fLBA * Lr, where JlBA ~ {ALBA/AArp220/Arp299A) X 
{LlR,LBA/Liji^Arp220/Arp299A) with logLifl,Arp220 = 12.27 

and \ogLiB^Arp299A = 11.53. fLBA was typically in the range 
1 — 2, and was never larger than a factor of 5, indicating 
that the comparison sample is well-matched to the LBAs in 
terms of their SFRs and size of the physical regions probed. 
Next, we calculate the ratio of observed to expected radio 
luminosity for both LBAs. 

For J0150-I-1308, the radio luminosity observed is, re- 
spectively, 1.7x and 2.8x higher than that expected based 
on Arp220 and Arp299A. For J092H-4509, and assuming 
that the source is resolved (see Section 3), the observed lu- 
minosity is, respectively, 2.9 x and 4.8 x higher than that 
expected based on Arp220 and Arp299A. If we assume that 
J0921-I-4509 is unresolved, its smaller area leads to observed 
luminosities that are higher than expected by factors of 4.7 
and 7.7. Similarly, if we compare with the brightest ra- 
dio nucleu s in Mrk273, Nl with L1.7 ghz « 7 x 10^** 
1 



erg 



s " Hz~^ (jBondi et al.ll2005^ . the observed luminosities are 
again higher by factors of 2-3 for both LBAs. To summa- 
rize, our comparison shows that the radio luminosities of the 
LBAs are always higher than expected by factor of 2-3 for 
J0150-hl308 and a factor of 3-8 for J0921-h4509 suggesting 
that more than just star formation may be contributing to 
the radio emission observed. 



4.3 The starburst interpretation 

In lOverzier et al.l ((20091) we argued that one explanation for 
the position in the BPT diagram of composite LBAs such as 
J0921-I-4509, together with their relatively low Lhci/Lfir, 
was that these are post-starbursts where the O stars were 
deficient and ionization by shocks (rather than photoioniza- 
tion by O stars) was very significant (our so-called "aged- 
starburst" hypothesis, §4.2.1 in that paper). However, con- 
firmation of shock heating in the form of enhanced [SII] 
and [OI] line strengths was not found. New COS data on 
J0921-I-4509 i mply a young age of ~3 Myr for the cen- 
tral starburst (jHeckman et al.l (|201ll ) and Borthakur et al., 
in prep.), ruling out that model. These are young star- 
bursts where the ionizing luminosity is at least an order- 
of-magnitude larger than the rate of kinetic energy supplied 
by supernovae. While t his does not mean tha t shocks may 
not be contributing (see lHeckman et al.l (|201in for evidence 
of high velocity outflows), there is certainly no evidence for a 
weakness or an absence of ionizing radiation from stars that 
would lead to shock-dominated emission lines if this were 
an aged starburst. The starburst interpretation also fails to 
explain the excess X-ray emission observed in this source. 



^ It should be noted that lDownes &: EckartI l|2003) found evidence 
for a weak AGN in the West nucleus of Arp220. If this is the case 
then the radio excess we detect would either be the result of 
many additional RNSe/SNRs, though that would be unlikely, or 
the result of a stronger AGN. 
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unless an exotic sta rburst contribu tion to the hard X-ray 
emission is invoked (jjia et aLll2011f) . 



4.4 The AGN interpretation 

A more likely interpretation is that (some) composite LBAs 
host an AGN, consistent with the large excess compact 
radio emission observed on VLBI scales (in the case of 
J0921-I-4509). Unfortunately, we do not have sufficient reso- 
lution and sensitivity to look for further confirming evidence 
by resolving the core-jet structure in the radio, and the non- 
detections at 5 GHz imply that these low-luminosity sources 
may have signifi cantly steeper spectra than found in radio- 
loud AGN (e.g. IPrandoni et al.ll2010l : iGiroletti fc Panessal 
[2OO9I). On the other hand, there are now four independent 
indications that an AGN may be present: (1) the composite 
position in the BPT diagram, (2) the high VLBI luminosity, 
(3) the c ompact radio s ize, and (4) the hard X-ray excess 
found bv lJia et al.l (|201ll . see also Fig. [3|. The offset into the 
AGN/composite region of the BPT diagram could then be 
the result of the strong ionizing radiation from the starburst 
with a small contribution from an obscured AGN (as well 
as some contribution from shocks). 

If we assume that the VLBI detection in J092H-4509 is 
due to a radio AGN, we can compare it to typical radio-quiet 
quasars and Seyfert ga l axies . Taking Lx,2-io ~ 2.3 x 10*^ 
erg s~^ from ljia et al.l (|201ll ) and applying a factor of 2.86 
to convert to Lx, 0.2-20 we find Lr/Lx ~ 1.3 x 10^*, about 
10 X higher th a n infe rred from the Lr/Lx relation from 
iLaor fc Beharl (J2008l ) who find < Lr/Lx >~ 10"^ (with 
a la scatter of about a factor 3 in Lr) for radio-quiet 
quasars. We com pare this sourc e to th e BH "fundamental 
plane" (FP) from iMerloni et all (J200l M03). Even though 
J092 1-1-4509 was not detected at 5 GHz, we have extrapo- 
lated our 1.7 GHz detection to 5 GHz assuming a spectral 
index of -0.1 (consistent with the (5(t) upper limit for this 
source in the 5 GHz observation). At face value, the result 
suggests a BH mass of a few times 10^ Mq . Much higher val- 
ues are unlikely given the constraint on the radio luminosity, 
although there is significant uncertainty in the BH FP de- 
pending on the source of t he radio and X-ray emission and 
the ac cretion state (e.g. s ee Plotkin et al.ll201ir . In a recent 
paper. Ide Gasperin et al.l (|201l|) find that a sample of low lu- 
m inosity Type II AGN and LINERS do not lie on the BH FP 
of lMerloni et al.l (J2003l ). Instead, based on their BH masses 
inferred from stellar velocity dispersions, the M03 relation 
would overpredict the Mrh of their so urces by up to two or- 
ders of magnitude. In ljia et al.l ()201ll ) it was estimated that 
the candidate Type 2 AGN in the composite LBAs may have 
BH masses in the range of 10^ — 10® Mq if these objects are 
radiating at the Eddington limit, but they c ould be higher 
for su b-Eddington BHs. On the other hand, lOverzier et al.l 
((20091) showed that the central mass concentration on a scale 
of ~100 pc is consistent with a "proto-bulge" of order lO'' 
Mq in stars. Assuming a standard bulge-black hole mass 
scaling relation, it is hard to imagine the presence of a black 
hole substantially more massive than ~ lO'' Mq. The pres- 
ence of a compact radio core has at least made it convinc- 
ingly clear that an AGN is present in J0921-I-4509, but it is 
impossible to give better constraints on the BH mass than 
lO^-lO'^ Mq. 

We can compare the results from our small sample with 



the results of IParra et al.l (|20ld ) who analyzed the radio 
AGN fraction among a much larger sample of relatively 
nearby, "warm" IR-luminous galaxies (wLIRGs) across the 
BPT diagram, also using VLBI. In some respects, their sam- 
ple can be considered to be the m ore obscured equival ent of 
our UV-selected LB A sample fsee lOverzier et al.ll2009l ). The 
wLIRGs and LBAs are shown in Figure[2] (small dots), where 
sources with large circles indicate VLBI detections. While 
most of the wLIRG sources with optical spectra indicative 
of a dominant AGN were detected in VLBI, most (~ 80%) of 
the sources in the AGN-starburst composite region of the di- 
agram were, in fact, not (see also lCorbett et al.ll2003l '). This 
situation is similar to the non-detections for three of our four 
composite LBAs. One may naively expect that the chance of 
detecting an AGN increases significantly as one moves from 
pure Hll-like objects to pure AGN-like objects across the 
BPT diagram, contrary to what is observed. One explana- 
tion for the lack of VLBI detections among the composite 
LBAs and wLIRGs could be that the high gas columns of 
the dense nuclear gas associated with the starbursts leads to 
significant free-free absorption, which would especially affect 
these composite systems. In this case, it would be better to 
observe at higher frequencies in order to offset the A^ de- 
pendence of the free-free opacity, although we note that we 
did not detect anything at 5.0 GHz either (including the 
known 1.7 GHz source J0921-f4509). On the other hand, 
even with a relatively high Lr/Lx of 10~* as found for 
J0921-I-4509, the remaining composite LBAs sources would 
have been hard to detect given their hard X-ray luminosities 
(about a factor of 10 lower). Not only are the radio AGN 
expected to be intrinsically faint, the relatively high redshift 
of our sources {z ~ 0.2) complicates their detectability with 
VLBI. Altern atively, given AG N variability as observed in 
some LLAGN (jBell et al.ll201lh . or the episodic accretion ex- 
pected in some cosmologic al galaxy formation models (e.g. 
iHopkins fc Quataertll20ld ) it is perhaps not surprising that 
our very small sample observed produced only one clear de- 
tection of an accreting BH in the radio as the X-ray obser- 
vations were performed at a different epoch. 



5 SUMMARY AND FUTURE PROSPECTS 

This paper presents the results of a pilot program to observe 
LBAs with composite spectra intermediate between pure 
starbursts and Type II AGN for comparison with objects 
whose optical emission-line ratios did not appear to suggest 
the presence of an obscured AGN. The VLBI data obtained 
with the EVN were combined with data previously obtained 
in the X-ray from XMM-Newton, optical from SDSS and 
mid- and far- IR from Spitzer to further constrain the na- 
ture of these objects that are of special interest given their 
demonstrated similarities to galaxies that were common only 
at high redshift. While it might be expected that compos- 
ite objects are detected more frequently than pure Hll-like 
regions only ~ 20% of our composite LBAs were detected, 
possibly due to free-free absorption. 

Three sources, J0150-H308, J0921-^4509, and 
J1029-I-4829 were detected above 5(t, and all had brightness 
temperatures close to or in excess of the typical upper-limit 
for thermal emission from star-forming regions. In addition, 
J0150-(-1308 and J0921+4509 had radio luminosities 2-8 
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times above that expected from multiple radio supernovae 
and supernova remnants in the starbursts Arp220, Arp229A 
and Mrk273 after a careful normalization based on the 
different resolutions and total SFRs. 

In particular, the composite LBA J0921+4509 displayed 
a compact core, most likely unresolved as confirmed by 
Monte Carlo simulation, emitting around 10% of observed 
VLA flux. The source also possessed an excess Lr expected 
from the Lh/Lx ratio of radio-quiet AGN, while its Lx was 
found to be high rela tive to the value expected based on its 
SFR (|jia et al.ll201ll ). This suggests a source of radio and 
X-ray emission beyond that expected of a pure starburst. 
We conclude therefore that this excess is most probably due 
to an obscured AGN with an estimated BH mass of ~ 10^~^ 

Mq. 

As stated in §1, one of our main motivations for probing 
BH formation in these LBAs as local analogs of LBGs is the 
question of whether the young starburst-dominated galaxies 
seen at high redshift already carry the "seeds" of supermas- 
sive BHs detected in the local universe. If so, these small-to- 
moderate mass BHs at high redshift may reveal their pres- 
ence through AGN activity. It is currently not clear whether 
BHs would be able to accrete very efhciently in these young 
and dense starburst dominated regions, or whether a delay 
should exist between the onset of a ce ntral starburst and 
the accretion onto the c entral BH (e.g. iNorman &: Scovillj 
ll988l : lDavies et al]|2007l '). Taken together, our X-ray and ra- 
dio results shown here indicate that LLAGN can co-exist 
in young, dense starbursts, but our sample is currently too 
small to put strong constraints on this timing argument. 

It is important to point out that at radio wavelengths, 
the contribution of these AGN is still only a fraction of the 
total radio flux density, making the identification of such 
sources based on radio data practically impossible much be- 
yond z ~ 0.2 (see below). Even in the X-rays, the com- 
bined luminosit y from SFR and candidate AGN in LBAs 
(|jia et al.ll201ll . and Fig. 3) is lower or comparable to that 
detected in t he deepest X-ray stacks of star-formi ng galaxies 
at z ~ 3 - 4 iLaird et al.ll2006l : IZheng et al.ll2010l '). illustrat- 
ing that individual detections are likely out of reach for cur- 
rent X-ray observatories. Also discussed in this paper, the 
interpretation of optical emission line diagnostic diagrams 
for identifying weak AGN is not trivial either, especially in 
the case of shocks or due to the limited phy sical resolution 
achieved in high redshift data (see Fig. 12 in lOverzier et al.l 
I2OO9I . for a large compilation of BPT information from low 
and high redshift galaxies). These limitations indicate the 
important role that nearby galaxies still play in studying 
the AGN-starburst (or BH-host galaxy) connection. 

Our pilot program of searching for AGN in the most 
extreme UV-selected starbursts in the relatively nearby uni- 
verse demonstrates the unique power of VLBI to assess the 
prevalence of accreting BHs and SN complexes in young, 
star-forming galaxies. With this goal in mind, we conclude 
with investigating the prospects for future VLBI detec- 
tion experiments. In Fig. |4] we show the predicted 1.4 GHz 
flux densities for (weak) radio AGN as a function of red- 
shift up to 2; = 1. The blue, green and red shaded ar- 
eas show the 1.4 GHz flux densities expected for a com- 
pact core having Lr/Lx ratios of 10~^ (i.e. typical radio- 
quiet AGN), 10^^ (i.e. that measured for the LBA composite 
J092H-4509), and 10^'^. The upper boundary of each region 



1000 



"3. 




0.4 0.6 
Redshift 



Figure 4. Predicted flux densities for low luminosity radio AGN 
as a function of redshift. The blue, green and red shaded areas 
show the 1.4 GHz flux densities expected for a compact core hav- 
ing Lji/Lx ratios of 10^^ (i.e. typical radio-quiet AGN), 10~* 
(i.e. lOx more luminous), and 10""^ (i.e. lOOx more luminous). 
The upper boundary of each region assumes Lx ~ 5 X 10''^ erg 
s~^ (i.e. comparable to that of J0921-I-4509), while the lower 
boundaries arc for sources that arc a factor 10 less luminous in the 
X-rays (i.e. comparable to the other composite LBAs studied in 
this paper). Our detection of the probable AGN in J0921-I-4509 is 
indicated by the star, while two undetected composite LBAs are 
shown as upper limits. Horizontal lines indicate typical EVN 5o" 
sensitivities achieved in a 1 hr (dotted line) and 4 hr observation 
(solid line) both at 1 Gbps recording, as well as the projected 4h 
sensitivity of the EVN with new telescopes becoming available in 
the next few years (dashed line). 



assumes Lx ~ 5 x 10'''^ erg s"'^ (i.e. comparable to that of 
J0921-I-4509), while the lower boundaries are for sources that 
are a factor 10 less luminous in the X-ray (i.e. comparable 
to the other composite LBAs studied in this paper). Hor- 
izontal lines indicate typical EVN 5cr sensitivities achieved 
in a 1 hr (dotted line) and 4 hr observation (solid line) both 
at 1 Gbps recording, as well as the projected 4 hr sensitivity 
expected in the next few years (dashed line) . A source simi- 
lar to J0921-I-4509 both in Lx and Lr/Lx can be detected 
to a redshift of z ~ 0.4, while it could be seen out to 2 ~ 0.8 
if its radio-loudness {Lr/Lx) was lOx higher. In contrast, 
the prospects for sources that are radio-quiet and a factor 
ten less luminous in the X-ray are grim: such sources could 
potentially be detected only out to 2 ~ 0.1. 

The recently increased sky coverage in the UV 
(GALEX), optical (SDSS) and IR (WISE) wiU allow the ef- 
flcient selection of exotic starbursts, composites, and AGN 
at relatively low redshifts suitable for VLBI. A systematic 
survey of their pc-scale radio emission offers a unique chance 
of gaining new insights into the complex nuclear processes 
occurring in galaxies over a very wide redshift range. 
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